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Abstract. Given a centred distribution, can one find a time- homogeneous 
martingale difTusion starting at zero which has the given law at time 1? 
We answer the question affirmatively if generalized diffusions are al- 
lowed. 



1. Introduction 

Consider a distribution with finite mean on the real hne. Can this dis- 
tribution be recovered as the distribution of a diffusion at time 1? Clearly 
not if the support of the distribution is disconnected since diffusions are 
continuous, but the answer is yes, and in many ways, if the target distri- 
bution is sufficiently regular. What if we also require our process to be a 
martingale or to be time-homogeneous? Again, easy constructions exist (see 
Remark 12.41 below) to show that a suitable process exists. But, what if we 
require our process to be both time-homogeneous and a martingale? The 
original motivation for us to study this problem comes from a calibration 
problem in mathematical finance. 

Our approach involves the speed measure of a diffusion and time-changes 
of Brownian motion, and makes no assumptions on the regularity of the 
target law. Indeed, to allow for target distributions with arbitrary support 
the natural class of processes to consider is the class of generalized diffusions 
(sometimes referred to as gap diffusions) discussed below. Our main result 
is to show that within this class of processes there is a time-homogeneous 
martingale with the given distribution at time 1. 

In Section [2] we introduce the class of generalised diffusions as time- 
changes of a Brownian motion. We also formulate our main result. The- 
orem 12. 3|, which states that given a distribution on M with finite mean, 
there exists a generalised diffusion that is a martingale and that has this 
distribution at time 1. In Section [3] we collect some general results about 
time-changes and generalised diffusions. In Section H] we study a discrete 
version of the inverse problem. For a given distribution with mass only in a 
finite number of points, we show that there exists a time homogeneous mar- 
tingale Markov chain with the given distribution at time 1. In Section [5] we 
consider the case of a general distribution on the real axis. By approximat- 
ing the distribution with finitely supported distributions, the existence of a 
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solution to the inverse problem is obtained, thereby proving Theorem 12 .31 In 
Section [H we apply our results to solve a calibration problem in mathemat- 
ical finance. In Section [7] we characterize the speed measures for which the 
corresponding generalised diffusion is a (local) martingale. Finally, Section [8] 
concludes with some open problems. 

2. Construction of generalised diffusions 

In this section we construct, following [6j (see also [llj and [12]), time- 
homogeneous generalised diffusion processes as time-changes of Brownian 
motion. The time-change is specified in terms of the so-called speed measure. 

Let be a nonnegative Borel measure on the real line; v may be finite 
or infinite. Let Bt be a Brownian motion starting at xq, and let be its 
local time at the point x up to time u. Recall that (a.s.) is continuous 
in (n, x) S [0, oo) x M and increasing in u. Define the increasing process 

r„ := / Llu{dx), (2.1) 
JR 

noting that G [0, oo], and let 

At := inf{n : F„ > t} (2.2) 
be its right-continuous inverse. The process 

'■= 

will be called a generalised diffusion with speed measure v. 

Example 2.1. If u{dx) = Jf^ for some continuous non-vanishing function 
a, then Xt is a weak solution of 

dXt = a{Xt)dWt. 

In this case, Xt is a diffusion. This example is the motivation for calling the 
measure 'speed measure', but note that i/ rather measures the inverse of 
the speed. 

Remarks 2.2. (i) Almost surely, oo for every x as u — )• oo; hence 

if v is non-zero, then F^ — ?> oo a.s. as n ^ oo, and thus At < oo and 
Xt is well-defined for every t G [0, oo) a.s. However, we have to 
exclude the exceptional case = 0, when F^ = for every u and 
At = oo for every t > 0, so Xt is not defined. (For technical reasons, 
we allow = when discussing At, but we always assume v ^ 
when considering Xt, sometimes without explicitly saying so.) 
(ii) F„ is left-continuous (by monotone convergence) and continuous at 
every u such that F^^ < oo (by dominated convergence); hence 
li 1-^ F„ is continuous everywhere except that there might exist a 
single (random) finite uq where F^ jumps to +oo: Tuq < oo but 
F^t = oo for all u > uq. (For example, this happens if v is an infinite 
point mass, see Example 13.61 below.) 
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(iii) By ()2.2p and the left-continuity of r„, for all t, u > 0, 

At<u <^ Fu > t. (2.3) 
Equivalently, u < At <J=^ t > Tu. It follows that < t < TAt+ 



for all t, and thus, by (ii) = t for all t such that TAt+ < oo; see 
further Lemma 17.41 
(iv) Tu is finite for all u a.s. if and only if is a locally finite measure 
(also called a Radon measure), i.e., iy{K) < oo for every compact K. 



In this case, Tu is continuous by (ii) , and is a continuous additive 



functional ( CAF) of B; conversely, every CAF of Brownian motion 
is of this type, see [TOl Theorem 22.25] and [HI Chapter X]. 

(v) Although Bq = xq, Aq may be strictly positive, and in general Xq ^ 
xo, see Lemma [3^ below for a precise result. In j6], Xq is defined to 
be xo, but this has the disadvantage of making Xt possibly not right- 
continuous at t = 0. We follow here instead the standard practise 
of considering right-continuous processes. We define Aq- = and 
Xq- = -Bao- = and thus allow the possibility that Xg- ^ Xq. 

(vi) We let {Tt)t>o denote the standard completed Brownian filtration. 
Then each At is a (J"j)-stopping time, and Xt is adapted to the 
filtration {Qt) = {^At)-, i ^ 0—. (In particular, t/o- = -^o is trivial.) 

In the sequel, we let "stopping time" mean (J-t)-stopping time 
unless we say otherwise. 

(vii) Even though the process Xt is constructed as a time change of Brow- 
nian motion, it is in general not necessarily a martingale or even a 
local martingale; see Section[7]for a detailed discussion. However, we 
are mainly interested in cases in which Xt is a local martingale, and 
preferably an integrable martingale. Recall the convention just made 
that Xq- = xq while Xq may be different. We say that {Xt)t>o- is a 
(local) martingale if {Xt)t>o is a (local) martingale (for the filtration 
(Gt)) and, further, KXq = Xq^. (This is equivalent to the standard 
definition interpreted for the index set {0—} U [0, oo).) 

Our main result is that any given distribution with finite mean can be 
obtained as the distribution of Xi for some such generalised diffusion with 
a suitable choice of speed measure i^. 

Theorem 2.3. Let fi be a distribution on the real axis with finite mean Jl = 
f^xiJ,{dx). Then there exists a generalised diffusion X such that Xq- = Jl, 
{Xt)o^<t<i is a martingale, and the distribution of Xi is /i. Furthermore: 

(i) Xq =]! if and only ifJiG supp /i. 

(ii) lKAi= Var(/i) := J^{x — Jl)'^ fi{dx) < oo. In particular, K Ai < oo if 
and only i/Var(/u) < oo. 

It follows from Theorem 17.91 that, actually, (Xj)o-<t<oo is a martingale. 

Remark 2.4. Theorem 12.31 guarantees the existence of a time-homogeneous 
(generalised) diffusion which is a martingale and has a certain distribution at 



4 



ERIK EKSTROM, DAVID HOBSON, SVANTE JANSON, AND JOHAN TYSK 



time 1. Note that the problem is much easier if some of these requirements 
are dropped. For example, for a given distribution ^, one can find a non- 
decreasing function / : M ^ M such that has law /i (in particular, 
/ = o where is the cumulative distribution function associated 
with n and $ is the normal CDF) and then the process Xt = f{Bt) is time- 
homogeneous with distribution /i at time 1. Further, if has a density which 
is strictly positive and differentiable, then X is a time-homogeneous diffusion 
with dXt = a{Xt)dBt + h{Xt)dt where a := f o /"^ and 2b := f" o f-\ 
Typically, however, it is not a martingale. 

Similarly, with / as above, the process Mt = E(/(i?i) | Tt) is a martingale 
with distribution /i at time 1, but typically it is not time-homogeneous. 
Again, in the regular case, we have Mt = g{Bt,t) for some function g solving 
the heat equation, and dMt = g'{g^^{Mt,t),t) dBt is a time-inhomogeneous 
martingale diffusion, where g~^{-, t) is the space-inverse of g{-,t), and g'{-,t) 
denotes the spatial derivative. 

Remark 2.5. In [9j, a similar inverse problem is studied. Instead of allowing 
for generalised diffusions, the authors consider the case of diffusions with 
regular diffusion coefficients obtaining an approximate solution to the inverse 
problem using variational techniques. In [3] the authors solve a related 
inverse problem in which the goal is to construct a generalised diffusion Y 
such that y is a martingale and such that the distribution of 1^^ is /x, where 
te is an independent random exponential time. Stopping at an independent 
exponential time is more tractable than stopping at a fixed time. 

Remark 2.6. One interpretation of our results is that we construct a stop- 
ping time T = Ai such that B^ = B^^ = Xi has law fi. The general problem 
of finding a stopping time T such that Bt has a given distribution is known 
as the Skorokhod stopping problem, see e.g. |13l Section 5.3], where other 
constructions of such stopping times are given. 

Example 2.7. As a very simple case, if fj, is the normal distribution N{xq, cj^), 
we may take di' = a~'^dx, a constant multiple of the Lebesgue measure. 
Then F^ = a~'^u and At = a'^t (both non-random), so Xt = B^2t, cf. Ex- 
ample [2Tj In general, however, it seems difficult to find i' explicitly. 

3. Preliminaries 
Recall that the support of a measure on M is 

supp V := is open and u{U) = O}. 

In other words, x G suppi/ if and only if every neighbourhood of x has 
positive measure. 

Similarly, we define the infinity set suppoo by 

suppoo z^ := M \ C M : [/ is open and u{U) < ooj. 
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In other words, x £ suppoo if and only if every neighbourhood of x has 
infinite measure. Thus, suppoo = if and only if v is locally finite. By 
definition, supp v and suppoo are closed subsets of M. 

If S is any Borel set in M, we let Hs := inf{t > : Bt G S} denote 
the hitting time of S (for the Brownian motion Bt). Note that if 7^ 0, 
then Hs < 00 a.s. We will only consider cases when S is closed, and then 
Bhs G S. For X € M, we write Hx for H^^}- 

Lemma 3.1. If v ^{), then a.s. Xf = B^t G suppj^ for all t >0. 

Proof. If s > and Bs ^ suppz^, then there exists e > such that Bu ^ 
suppz^ for n G [s — e, s + e], and then is constant for n G [s — e, s + e] for 
each X G supp i^; hence is constant for u G [s — e, s + e], and thus At s 
for all t > 0. 

Similarly, if s = and Bg ^ supp v, then there exists e > such that 
F„ = for u G [0, e]. Therefore At > = s for all t>0. □ 

Lemma 3.2. IfT is any finite stopping time for Bt, then a.s., for all e > 0, 

Consequently, there is for every e > a.s. a (random) 5 > and an open 
set U containing Bt such that L^+e ~ — ^ /'^'^ '^^^ x (z U. 

Proof. The first claim is an immediate consequence of the strong Markov 
property and the fact that > a.s. for a Brownian motion started at x 
and any e > 0. 

The second claim follows since L is continuous. □ 

Lemma 3.3. Aq = Hgy^ppu, the hitting time of the support of v, a.s. In 
particular, if xq G suppj^, then Aq = Q and = xq a.s., hut if xq ^ suppz/, 
then ylo > and Xq ^ xq a.s. 

Proof. If u < H = Hsuppp, then F„ = 0. On the other hand, by Lemma [3. 21 
for any e > 0, L^"^^ > a.s. and, moreover, L^^+e > for x in a neighbour- 
hood U of Bh- Since Bh G suppz^, we have ^{U) > and thus Tn+e > 0. 
In other words, a.s. F^ > for all u > H. 

The definition (12. 2p of At now shows that Aq = H a.s., and the result 
follows, recalling Lemma l3.1[ □ 

Lemma 3.4. We have KXq = xq = Xq^ if and only if neither suppz^ C 
(— oo,xo) nor supp C (xo, 00). 

Proof. If supp J/ intersects both [xo,oo) and (— oo,xo], then BuAH^uppi, is a 
bounded martingale and thus EXq = ^Bn^uppi, = Bo = xq. 

On the other hand, if suppi^ C (— oo,xo), then Lemma [3T] implies Xq < 
Xq a.s., so KXq < Xq. Similarly, if suppz^ C (xo,oo), then EXq > xq. □ 
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We define, for given v and xq, 

x+ := inf{x > xq : x G suppz^}, (3-1) 

x_ := sup{x < xq : X S suppi^}, (3-2) 

:= inf{x > xq : x G suppoo J^}, (3-3) 

x'^ := sup{x < Xq : X G suppoo 1^}- (3-4) 



Note that these may be ±00 (when the corresponding sets are empty). In 
general, 

—00 < x'^ < x_ < xo < x+ < x'^ < 00. 
It fohows from Lemmas 13.31 and 13.41 that we have the fohowing cases: 

(i) If Xq G suppi/ (i.e., x_ = x+ = xq), then Xq = xq. 

(ii) If suppi/ C (xo,oo) (i.e., x_ = —00 and xq < x+ < 00), then 
Xq = x+. 

(iii) If suppz^ C (— oo,xo) (i.e., x+ = +00 and —00 < x_ < xq), then 

Xq = X-. 

(iv) Otherwise (i.e., if —00 < x_ < xq and xq < x_|_ < 00), then Xq G 
{x_,x+}, with the unique distribution satisfying EXq = xq. 

Lemma 3.5. Let H = -ffsuppoo f hitting time for Bt of the infinity set 

of V. Then a.s. L^ = 00 for all u> H , and thus At < H for all t >0. 
On the other hand, L^ < cxd a.s. for all u < H. 

Proof. If i7 < 00, then Bh G suppoo i^- For any e > 0, by Lemma [3.21 
L|^^^ > 6 > for X in a neighbourhood U of Bh and some 6 > 0. Since 
Bh G suppoo z^, we have viU) = 00 and thus Tn+e > Ju Lf[^^i'{dx) = 00. 
The definition (12. 2p of At now shows that At < H a.s. 

If u < H, then K := {Bg : < s < u} is a compact interval disjoint from 
suppoo i^, and thus I'^K) < 00. Since x 1— > is continuous and vanishes 
outside K, we have L^ = jj^ L^v{dx) < 00. □ 

Example 3.6. Let be an infinite point mass at xi G ffi. Then suppz^ = 
suppoo = {xi}. If is the hitting time of xi, then T^ = for u < f^xi, 
cf. Lemma 13.31 and its proof, but a.s. F^j = 00 for u > H^^ by Lemma 13.51 
Hence, a.s.. At = Hx^ and Xt = xi for every t > 0. 

In particular, if xi = xq, then Hx^ = a.s., and thus At = Q a.s. for all 
t > 0. 

More generally, if z/ is a measure such that i^(5) = or I'^S) = 00 for 
every Borel set S, then suppoo = suppz^, and Lemmas 13.51 and 13.31 implv 
that a.s. At = Hguppu for all t > 0, so Xt = Xq is the first point of suppz/ 
hit by Bu. 

Lemma 3.7. Suppose that suppoo 1^ H (—00, xq] and suppoo H [xq, 00) are 
both non-empty. Then BaiAu £ [x!^, x^^] for all t,u > 0, and thus (-B^t am)m>o 
is a bounded martingale for each fixed t. 
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Proof. Recall from ()3.3|) ~ (I3.4|) . By assumption, — oo < x'^ < xo < x'^ < 

Then < 5„ < for all u < H; 
thus Lemma 13.51 implies that x"^ < B^AAt ^ 2:+ for any n > and t > 0. 
Finally, {BAtAu)u>o is a martingale since At is a stopping time. □ 

Lemma 3.8. IfEAt < 00, then EXt = xq and Yar Xt = E{Xt - xq)^ = 
EAt < 00. 

Proof. This is an instance of Wald's lemmas, see e.g. |13l Theorems 2.44 and 
2.48]. □ 

Lemma 3.9. IfEAt^ < 00 for some to < 00, then (Xt)o-<t<to is a square 
integrable martingale. 

Proof. Since each A n is a bounded stopping time, if 0— < s < t, then 

E{BAfAn I J^As) = BAtAnAAs = Ba^Au (3.5) 

a.s. (see e.g. [TOl Theorem 7.29]), so BAtAn, t > 0—, is a martingale. Fur- 
thermore, by Wald's lemma, i.e., since {Bt — xq)^ — t is a martingale, 

E{BA,An-xof = E{AtAn). (3.6) 

It follows that for any fixed t < t^, E{BAtAn — xq)'^ <EAt< 00; hence the 
variables BAtAn, n > 1, are uniformly integrable, and thus BAtAn — > BAt 
in L^. If 0— < s < t < we thus obtain, by letting n — > 00 in (j3.5p . 
E{BAt I J'as) = Ba^ a.s. 

Thus Xt = B^t, 0— < t < to, is an integrable martingale; it is square 
integrable by Lemma |3. 81 □ 

Note that the converse to Lemma 13.81 does not always hold: we may have 
Yar Xt < 00 also when EAt = 00. For example, this happens in Example l3.6l 
if xi 7^ Xq. We give a simple sufficient condition for EAt < c«. 

Lemma 3.10. Suppose i/iai suppoo i^H (— cxd, xq] and suppoo i^H [xq, cxd) both 
are non-empty. Then {Xt)t>o-! is 0, hounded martingale with EXt = xq and 
EAt = E{Xt - xof < 00 for every t>0. 

Proof. Lemma 13.71 shows that x"^ < Xt = BAt — ^+ ! fo'^ t >0, so (Xt) 
is uniformly bounded. 

For each n, A n is a bounded stopping time and (j3.6p holds. Letting 
n — > 00, we find E{Xt — xq)^ = EAt by dominated and monotone con- 
vergence, and thus EAt < max{(xo — x^)'^,{x'^ — xq)^} < 00. Finally, 



Lemma 13.91 shows that Xt is a martingale. □ 

We have defined At in ()2.2p so that it is right-continuous. The corre- 
sponding left-continuous process is 

At- := inf{u > : r„ > t}; (3.7) 



note that for t > 0, At- = lim.s/'t As, while Aq- = (as defined in Re- 
mark l2.^v)| above), and At- is a stopping time. It is possible that At- < At, 
i.e. that At jumps; this corresponds to time intervals where F^ is constant. 



8 



ERIK EKSTROM, DAVID HOBSON, SVANTE JANSON, AND JOHAN TYSK 



because moves in the complement of supp i^, so miless supp = M, it will 
a.s. happen for some t. However, the next lemma shows that there is a.s. 
no jump for a fixed t > 0. (Equivalently, for a fixed t > 0, there is a.s. at 
most one u > such that r„ = t.) 

Lemma 3.11. Let t he fixed with < t < oo. Then a.s. At- = At. 

Note that the result fails for t = 0, see Lemma 13.31 

Proof. Let e > 0. Since At- is a stopping time, Lemma 13.21 shows that 
a.s. there exists a neighbourhood U of Ba^- and some (5 > such that 
L\^_j^^ > L\^_ + 5 for X G [/. Further, since At-ij^ — > At- as n — )• oo. 
Lemma l3 . II implies that a.s. -Byit- £ suppi^ and thus ^{11) > 0; hence either 
^At- = oo or 

Ta^.+c > Ta,- + Siy{U) > Ta,. . (3.8) 

If TAt-+e < OO, then r„ is continuous at u = At-, see Remarks I2.^ni)[ 
and it follows from (13. 7p that r^^_ > t (and actually r^^_ = t); thus (j3.8p 
yields TAt-+e > t. This is trivially true also when TAt-+s = oo. 

Thus, a.s., TAt_+e > which implies that that At < At- +e. Since e > 
is arbitrary, and At- < At, the result follows. □ 

When considering several speed measures Un, we use n as a superscript 
to denote the corresponding T^, A^ and X"; we use always the same Bt. 

If 5 is a topological space (in our case M or an interval in M), we let Cc{S) 
denote the space of continuous functions S" — > M with compact support, and 
C^{S) the subset of such functions S — )• [0, oo). 

Lemma 3.12. Let u,ui,U2t ■ • be a sequence of measures on M. Assume 
either 

(i) /jj if dun — )■ J^(p du as n ^ oo for every (p G (R) , 
or, more generally, 

(ii) there exists an interval (a, b) with —oo <a<XQ<b<oo such that 
j ipdvn J ip dv as n oo for every ip G C^{a,b), and also for 
every ip G C^(M) such that p{a) > or p[b) > 0. 

Then, for each t > 0, Af — > At a.s., and thus, if v,vi,U2, ■ ■ ■ are non- 
zero, Xt" — Xt a.s., where X,X^,X'^, . . . are the corresponding generalised 
diffusions constructed from the same Brownian motion. 

Proof. The local time G C(^(M), as a function of x, for every u > 0. In 
(i) we thus have, for every > 0, 

Tl= j LlMdx) ^ [ Ly{dx) = r„. (3.9) 

In (ii), let H = H^g^ i,j be the hitting time of {a,b}. u < H, then the 
support of is contained in (a, 6), so LJ; G C^{a,b) and L" — )• as in 



CAN TIME-HOMOGENEOUS DIFFUSIONS PRODUCE ANY DISTRIBUTION? 9 



If u > H, then a.s. L^" > by Lemma [3]2j Since (p{x) := is continu- 
ous and Bh = a or Bh = b, the assumption shows that (13. 9p holds in this 
case too. 

Hence, in both (i) and (ii), F" for all u > except possibly when 

u = H. 

Let s = At- By (12. 2p . Tu> t for u > s and Tu < t for u < s. Further, if 
Tu = t for some u < s, then At- < u < s = At, which has probability by 
Lemma 13.111 Consequently, a.s. < t if u < s. 

Assume for simplicity that s = At < oo. (The case ^4^ = oo is similar.) 
If e > and s + e ^ H , then thus F"_,_^ F^+e > t. Hence, for sufficiently 
large n, r"^_^ > t and thus A^ < s + e. Since e > is almost arbitrary, 
it follows that limsup„^oo^" < s a.s. Similarly, considering s — e ^ H, it 
follows that lim inf„^oo ^j" ^ s a.s. Consequently, A^ At a.s. as n — )• oo. 

Finally, if At < oo, then X'^ = Ba^ B^t = Xt a.s. follows, since B^ is 
continuous. □ 

Lemma 3.13. Let a > b > xq and 6 > 0. Then there exists a constant 
C = C{a,b,xo,6) such that if¥{Xi >a)>5, then vlxQjb] < C. 

Proof. Assume for convenience xq = 0. By replacing a by (a + 6)/2, we may 
also assume that P(Xi > a) >6. 

Let H = Ha = inf{u : Bu = a}. By definition, Xi = Ba^, so if Xi > a, 
then H < Ai and thus F^ < 1. Consequently, F{Th < 1) > 5. 

The local time is a continuous function of x G M, and it is a.s. strictly 
positive on [0, a) by Ray's theorem [13, Thm 6.38] (a consequence of the Ray- 
Knight theorem which gives its distribution). Hence, Y = inf^.g[Q > 
a.s., and thus there is a constant c > such that ¥{Y < c) < 6. 

Hence, with positive probability Th < 1 and Y > c. However, then 

l>F//= / L%v{dx)> [ L'^H'^idx) >Yi^[0,b]>ci^[0,b], 
Jr Jo 

so iy[0,b] <l/c. □ 

Lemma 3.14. For every K > there exists k = k{K) > such that if u is 
a speed measure such that E < K , and further suppoo v n (— oo, xq] and 
suppoo V n [xo, oo) both are non-empty, then u[xq — 2K, xq + 2K] > k. 

Proof. We may assume that xq = 0. By Lemma 13. 7^ {BuaAi)u>o is a 
bounded, and thus uniformly integrable, martingale, closed by Bai = Xi. 
Let H = H^j-2K} be the hitting time of ±2K. Then 

P(F^ < 1) = FiH < A^) < P(sup \B^aaA > 2K) < ^ ^ = ^• 

(3.10) 

Let Y = max^. L^- • this is a finite random variable so there exists c > 
such that P(y > c) < 1/2. (Note that Y and c depend on K but not on 
u.) With positive probability we thus have both F^^^ > 1 and Y < c. This 
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implies, since L% = when Ixl > 2K, 

I-2K i'2K 

l<Tfj= j L}.i^{dx)< Yi^{dx) <cu[-2K,2K], 

J-2K J-2K 

and the result follows with k = c~^. □ 

4. The discrete case 

In this section we treat the inverse problem in a discrete setting. We fix 
points yo < yi < y2 < ■ ■ ■ < 2/n+i and consider discrete speed measures 
v = "^^=0 bi5y- , where 6a is a unit point mass at the point a and 6j takes 
values in [0,oo]. We assume that bo = bn+i = oo. We also fix a starting 
point xq G (yo,yn+i)- (We could for simplicity assume that xq = yi^ for 
some io G {1, 2, . . . , n}, but that is not necessary.) 

Given such a speed measure v and xq, we construct a generalised diffusion 
X as described in Section [2] above. By Lemma [3. 11 the process {Xt)t>o only 
takes values in the set {yi\^^Q ■ Moreover, since 6o = ^n+i = oo, the states 
yo and yn+i are absorbing, so X is bounded, and it follows from Lemma [321 
(or Theorem 17. 3p that X is a martingale; in particular EX^ = xq. Let 
Pi = P(Xi = yi) be the probability that X at time 1 is in state y^. Then 

< K < 1, Yll=oPi = 1> and we also have Yh^o ViPi = ^Xi = xq- 

This defines a mapping G from the set of speed measures above to the 
set of distributions with mean xq. More precisely, we write G(bi, . . . , bn) '■= 
(pO) • • • )Pn+i)) and note that G : i?" 11", where := [0, 00]*^ and 

n+l 71+1 

m := jvr = (vro, . . . ,7r„+i) G [0, 1]"+^ : ^ vr^ = 1 and ^ym = xq}. 

i=0 i=0 

Lemma 4.1. The function G : B"" —?■ 11"^ is continuous. 

Proof. This is a direct consequence of Lemma [3TT2J (Note that the possibility 
that one or several 6j = 00 is no problem when we verify condition (i).) □ 

We will use algebraic topology to show that G is surjective; see e.g. [U 
Chapter IV] for standard definitions and results used below. We begin by 
studying the sets B^ and 11*^. The set B^ is homeomorphic to the unit 
cube [0, 1]" with the boundary dB"' := [0, 00]" \ (0, 00)" corresponding to 
the boundary ^[0, 1]'' = [0, 1]" \ (0, 1)'^. We write 

n 

dB^ = [j (c/.oi?" U 5,ooi?") (4.1) 

where 5^,5" := {(6^)^ e B"" : bi = s}. 

The set 11" is the intersection An+i fl M^g of the simplex 

n+l 

A„+i := {(7r,)^+i : VTi > and ^ VTi = 1} 

1=0 
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and the hyperplane 

n+1 

in M"+2. Further, A„+i hes in the hyperplane L := { (vri)^^^ : vr^ = l} . 

Thus n" is a compact convex set in the n-dimensional plane L n Mx^ , which 
can be identified with R". Since yo < xq < yn+i, Mro contains an interior 
point of A„-|_i, so n" = A„+i n Mxq has a non-empty relative interior in 
L n Mxg ■ Thus n" is an n-dimensional compact convex set in L D Mxq and 
its boundary is 

n+1 

an" = n" n dAn+i = [j Sj-n", 

j=0 

where djU"" := G : ttj = O}. 

Consequently, both i?" and 11" are homeomorphic to compact convex 
sets in M" with non-empty interiors. Every such set is homeomorphic to 
the unit ball := {x G M" : \x\ < 1} via a homeomorphism mapping the 
boundary onto the boundary dD" = 5"^^. Thus there are homeomorphisms 
(5", (95") ^ (D",^"-!) and (n",9n") (L>",S'"-i). 

Lemma 4.2. The function G maps dB" into 511", ancf thus G : (-B", 9i?") — )• 

(n",9n"). 

Proof. If (ftj)" G djoB"^, then ^ suppz^, so Xi 7^ a.s. by Lemma [XT] 
and thus ttj = 0, so (vri)[}+^ G S^n". 

If G djooB"-, then G supp 00 i^. If further yj < xq, then Lemma [3^5] 
implies that Xi > yj a.s., and thus tTj = for < i < j and, e.g., (vrj)Q"'~^ G 
dj-iH". Similarly, if yj > xq, then tTj = for j < i < n -|- 1 and i'n'i)^'^^ G 

9j-+in". □ 

The homeomorphisms above induce isomorphisms of the relative homol- 
ogy groups i/„(S",aS") w i/„(n",9n") ^ Hn{D'^,S''-^) ^ 7L. The map- 
ping degree of the function G : (i3",5-B") — )• (11", 511") can thus be defined 
as the integer deg(G) such that the homomorphism : Hn{B'^,dB'^) — >■ 
i?„(n",5n") corresponds to multiplication by deg(G) on Z. More precisely, 
this defines the mapping degree up to sign; the sign depends on the orien- 
tation of the spaces, but we have no reason to care about the orientations 
so we ignore them and the sign of deg(G). 

Lemma 4.3. For any n>l, any yo < yi < ■ ■ ■ < y„+i, and xq G (yo, Vn+i), 
deg(G) = ±1. 

Proof. We use induction on the dimension n. We sometimes write G = Gn 
for clarity. For the induction step, we assume n > 2. The long exact 
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homology sequence yields the commutative diagram 

G* 

o = F„(n") — y /7„(n",an") H^-imn — > i^„_i(n") = o 

where the rows are exact; thus the connecting homomorphisms d are iso- 
morphisms, and the degree of G : (5", dB^) (11", dH"') equals the degree 
of the restriction G : dB^ — > 911". 

Assume xq < Dn- (Otherwise xq > yn > yi, and we may argue similarly 
using 9ioo-B" and ^qII".) We single out the faces dnooB"' and 9„+in" of 
the boundaries and define 0*5" := Ur=i ^iO^" U Ur=i^ ^ioo-B" and a*n" : = 
Ur=o^in". By the proof of Lemma [121 G : a„oo5" ^ 5„+in" and G : 

a*5" ^ a*n". 

We claim that the degree of G : dB^ — )• 511" equals the degree of 

G : (9„ooS",5noo5" n a*5") ^ (a„+in",a„+in" n a*n"). using home- 

omorphisms 95" ~ S""^^ and 911" ~ S*""^ that map the faces dnooB"" and 
9„+in" onto the upper hemisphere S"~^, this is an instance of the general 
fact that if F : S''^~~^ — )• S^~^ is continuous and maps S^~^ — > S^~^, then 
the degree of F : {Sl-\Sl-^ D ST^) (5!^-\5!^-^ n S'T'^) equals the 
degree of F : 5"-^ ^ S"~i. 

If (6i)" € 9„oo5", then has infinite point masses at both y„ and yn+i; 
with < yn < Un+i- By Lemma 13.51 we can ignore Vn+i and we obtain 
the same generalized diffusion X as with the speed measure X]"=o ^i^Vi ■ 
can thus identify dnooB^ with 5"~^ (based on the points (yj)o)- Further- 
more, there is an obvious identification 9„+in" = 11""^, and with these 
identifications, G : Snoo-B" — > 9„+in" corresponds to Gn-i ■ 5""^ 11""^. 
Moreover, the various boundaries correspond so that we have the commu- 
tative diagram 

5„_i(9,oo5",9,oo5"na*5") = 5„_i(5"-i,95"-i) 



Gn * 



where the rows are the isomorphisms given by these identifications. Hence 
the degree of Gn : (9„oo5", a„oo5" n 9*5") ^ (S^+iH", 9„+in" n 9*n") 
equals deg(G„,_i). 

Combining this with the equalities above, we see that deg(G„) = deg(G„_i) 
which completes the induction step. 

It remains to treat the initial case n = 1. In this case B^ and 11^ are 
intervals, and can be parametrized by bi and pi. It is easy to see that the 
mapping G : 6i i— > pi is strictly increasing, and thus a homeomorphism 
5^ — )• n^; hence G* : Hi{B^,dB^) — > 5i(n^,9n^) is an isomorphism, so 
deg(G) = ±1. Alternatively, we may use the first commutative diagram 
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above also in the case n = 1, replacing the homology groups Hn-i = Hq by 
the reduced homology groups Hq. The sets dB^ and 911^ contain exactly 
two elements each, and it is easy to see that G : dB^ dU^ is a bijection 
and thus : Ho{dB^) HQ{dIl^) is an isomorphism. □ 

Theorem 4.4. The function G is surjective, for any n > 1 and any < 
yi < ■■■ < Un+i and xq € (?/0) yn+i)- Consequently, the discrete inverse 
problem has a solution. 

Proof. An immediate consequence of Lemma l4.3| since a function {B^, dB"") — > 
(U",dW) (or, equivalently, ^ (Z)",^"-^)) that is not surjective 

has mapping degree 0. □ 

5. The general case 

In this section we study the inverse problem for arbitrary distributions 
on the real axis. To do this, we approximate the given distribution with 
a sequence of discrete distributions. For each discrete distribution we can 
find a discrete speed measure that solves the inverse problem according to 
Section m We then show that the sequence of discrete speed measures has a 
convergent subsequence, and that the limit solves the inverse problem. We 
begin with a lemma giving the approximation that we shall use. 

Lemma 5.1. Let fi be a probability measure on M with finite mean JI = 
J^x fi{dx). Then there exists a sequence fj,n, n>l, of probability measures 
with finite supports such that, as n ^ oo, 

(i) Hn ^ n weakly; 

(ii) inf supp fin inf supp fi; 

(iii) sup supp fin sup supp fi; 

(iv) each has the same mean Ji as fi. 

(v) /j. \x\ dfj.n{x) /]g \x\ dn{x). 

If further fi has finite variance Var(/i) = J^{x — Jl)'^ fj,{dx) , then fin can be 
chosen such that 

(vi) Var(;U„) ^ Var(;u). 

Proof. Let y be a random variable with distribution fi. First, truncate Y 
at ibn by defining 

:= (y An) V (-n). 
Then, letting all limits in this proof be for n — ?> oo, 

E IK - ^1 < IE(|y|; |y| > n) ^ 0. 
Next, discretize by defining 

Y::=-[nY;\. 
n 

Clearly, jF^' - < It follows that 

|iEy^'-Ey| <E[y^'-y| < i + E|y^ -y] ^ o. 
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Finally, we adjust the mean by defining 

F„:=y^'-E(y^'-y). (5.1) 

ThusEy„ = Ey = 71. 

Let fin = 'C(y„), the distribution of Yn- Then /i„ has finite support and 
(iv) holds by the construction. Furthermore, by (j5.ip . 

E|y„ -y| < 2E|y^'-y| ^0, (5.2) 

which implies y„ — ^ Y and thus (i). From (|5.2p we also have E |y„| — t- E |y|. 



which is (v) 



If inf supp/x = — 00, then (i) implies that inf supp^„ — >• — oo, so (ii) holds. 
Suppose now that inf supp = a> —oo. Ifn> |a|, then inf supp C{Y^) = 
a, and it follows that 

I inf supp /i„ -a\<- + lEy'' -Ey| ^ 0, 
n 

which shows that (ii) hold in this case too. 

The proof of (iii) is similar, mutatis mutandis. 

If /X has finite variance, thenEy^ < oo,andE|y^|2 = E{\Y\Anf ^Ey^. 
Taking square roots we find ||y^||2 11^ II 2- Minkowski's inequality yields 

|||>;^||2- lli^llbl < \\Yn-Y% < i + |E(y^'-y)| ^0. 

n 

Consequently, ||y„||2 ^ ||y||2, and thus Ey^^ ^ ^y^. Since Ey„ = Ey, 
this implies Var(l^) — )• Var(y), which shows (v). □ 

Proof of Theorem \2.3[ We may for simplicity assume that xo = Jl = 0. 

Let a_ = inf supp /X > — oo and = sup supp < +oo. Since /I = 0, we 
have a_ < < a+. Moreover, if a_ = or a+ = 0, then necessarily fi = 6o. 
In this case we may simply take u as an infinite point mass at 0; then At = 
and Xt = for all t > a.s., see Example 13.61 In the sequel we thus assume 
— oo < a_ < and < a+ < oo. 

Let Hn be a sequence of distributions satisfying (i) - (v) in Lemma 15.11 
The distributions have finite supports, and thus Theorem 14.41 shows that 
there exist speed measures Vn so that the corresponding generalised diffusion 
X" has distribution Hn at time 1. 

If < 6 < a+, choose a G (6, a+). Then a > 6 > and a < sup supp /i, 
so fi{a,oo) > 0. Since /x„ n, liminf„_!.oo /Xn(«> oo) > /x(o, oo), so for all 
large n, /i„,(a, oo) > ^^(a,oo) > 0. Lemma [3.131 applies and implies that 
z^n[0, 6] < C = C{b) for all large n, i.e., limsup„_^oo f„[0, 5] < oo for every 
6 < a+. 

By a symmetric argument, we also have limsup„_^oo i^nib, 0] < oo for every 
6 > a_. 

Choose sequences 6™ \ a_ and b'^ «+• On each interval [6™, 6™], the 
measures f„ are, as we just have shown, uniformly bounded if we exclude 
a finite number of small n. Since is compact, we may thus choose 

a subsequence of Un such that the restrictions to [b'^,b^] converge to some 
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measure Um- By a diagonal procedure, we can do this simultaneously for all 
m, which provides a subsequence of Vn such that (along the subsequence) 
z^n i'm on [6™, 6™] for every m. In the sequel we consider only this 
subsequence. 

It follows that A(/) = lim„^oo / di^n exists and is finite for every / G 
C+(a_,a+). Clearly, A is a positive linear functional on C+(a_., a^), so by 
the Riesz representation theorem there exists a Borel measure u on (a_, a+) 
with A/ = Jjg/d?. Thus f^fdi^n ~^ f f for all / € Cc(a_,a+). 

We define by adding infinite point masses at a_ and a+, if these are 
finite: = + oo • (5a_ + oo ■ Ja+ (where 6±oo = 0). 

If a_ is finite and / G C+(M) with /(a_) > 0, then f f du > /(a_)z^{a_} = 
oo. Furthermore, a" := infsupp^„ — )• a_, so f{a"L) > for all large n. The 
construction of i^n in Theorem 14.41 gives an infinite point mass at a"i, 
so / f di^n = C)0 for all large n. Thus / /df„ J f du = oo as n ^ oo. 
Similarly, J /df„ — > J / di^ = oo as n — > oo if /(a+) > 0. 

The assumptions of Lemma l3.12( ii) are satisfied, and thus a.s. A" — )• Ai 
and, if 7^ 0, X" Xi as n ^ oo. In particular, then Xf — > Xi in 
distribution, and since X^ has distribution /u„ and — ^ it follows that 
the distribution of Xi is /i. 

It remains to verify that the measure is non-zero. If a_ or a+ is finite, 
this is clear since by construction has a point mass there. 

If a± = ±oo, we use Lemma ISHT v) | which yields E \Xf\ = \x\ dfin{x) — ?> 

|x| dfi{x) < oo. Let K := sup„E \X^\ < oo. By Lemma [3.141 there exists 
K > such that i'n[—'^K,2K] > k for every n. Let / € Cj^( — oo,oo) with 
/ = 1 on [—2K,2K]. As shown above, J^fdi^n /r/^^ ~ luf^'^ ^ 
n — >■ oo. Since J^f dun > z^n[— 2K, 2Er] > k, this implies f diy > k > 0, 
and thus 7^ 0. 

This proves the existence of a non-zero speed measure z^ such that Xi has 
the desired distribution /i. We next prove that Xt is a martingale. 

We have shown that Xf Xi a.s. and, by Lemma I5.]](v)| E|X"| — ?> 
E[Xil. This implies E\Xf - Xi\ 0, i.e. X^ Xi in L^ee e.g. [13 
Proposition 4.12]). For each n, (-Ba5'au)u>o is by Lemma [3.71 a bounded 
martingale with limit Ba^ = A"", and thus Ba^^au = E(X" | J^u) for every u. 
As n — )■ 00, a.s. A" — t- Ai and thus Ba^au — ^ -Byi^/x^. Further, we have just 
shown — > Xi in L^, and this implies -Ba'j'au = E(Xf | J'u) E{Xi \ Fy) 
in L^. The two limit results both hold in probability, so the limits must 
coincide: Ba^au = E(Xi | Tu)- 

This proves that {BaiAu)u>q is a uniformly integrable martingale. Con- 
sequently, for any (/'t)-stopping time r, Ba^at = IE(Xi | Jv). In particular, 
for 0— < t < \, Xt = Bai = E(Ari I J-Ai), which proves that {Xt)t<i is a 



martingale. This completes the main part of the proof, and we turn to (i) 



and (ii) 



By Lemma [3^ ATq = /U if and only if G supp v. If /i ^ supp z^, then there 
exists an open set U with Jl £ U and z^(?7) = 0. By Lemma 13.11 Xi ^ U 
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a.s., SO JJ ^ supp/x. On the other hand, if /Z G suppi^ and U is any open 
set containing Ji, then I'iU) > 0. It is easy to see that there is a positive 
probabihty that Bt will remain in U until Fj > 1, and thus Xi G U. Hence 
/i(C/) = P(Xi G [/) > for any such U, so JI £ supp ^. 

If has finite variance, we may by Lemma IS.ll fv) assume that Var fin 
Var fi, and thus sup„ Var Hn < oo. Lemma 13.101 applies to every f„ and 
yields 

E^^ =E(Xi")2 = Var(/i„). 
Consequently, by Fatou's lemma, 

E^i = E lim A"^ < liminf E = liminf Var(/x„) = Var(/i) < oo. (5.3) 

Lemma 13.81 shows that KAi = Var(/i) in this case (so there is equality in 
()5.3p ). and also if Var(/i) = oo. □ 

Example 5.2. As an illustration of Theorem 12.3^ let Yt = where 

is a 3-dimensional Brownian motion with |-Bg^^| = Uq^ for some yo > 0. 
It is well-known that y is a local martingale bounded from below, hence 
a supermartingale, but it is not a true martingale, compare Theorem 17.91 
below. Moreover, Y can be represented as the solution of 

(5.4) 

for some standard Brownian motion W. Note that y is a diffusion by 
(|5.4p . and by Example 12.11 (together with a stopping argument), Y is the 
generalized diffusion with speed measure x~'^dx, x > 0. (Note that Yt tends 
to as i — > oo, but never reaches 0.) Being the reciprocal of a Brownian 
motion, Y has an explicit density, compare [21 Example 2.2.2]. The density 
of Yi decays like Cx~^ for large x, so only moments of order strictly less than 
three exist finitely. Moreover, the expected value of Yi is strictly smaller 
than the starting value yo- 

Theorem 12.31 provides the existence of a time-homogeneous generalised 
diffusion X which is a true martingale such that Xi and yi have the same 
distribution, and Xq = KYi < uq. Consequently, there are in this case two 
different speed measures that give rise to the same distribution at time 1. 
However, only one of the corresponding processes is a martingale, and they 
have different starting points. 

Example 5.3. A related example is when Y is the solution of 

(dYt = {l + Yt^)dWt, t>0, 
\Yo = 0. 

By Example 12.11 the diffusion Y is the generalized diffusion with speed 
measure {1 + x'^)~'^ dx . Again, y is a local martingale, but not a martingale, 
see Theorems 17.31 and 17.91 below. 
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Define h{x,t) = e g{x), where g is a smooth positive function satisfying 
g{x) = \x\ for l^l > 1 and g{x) > \x\ everywhere, and M is a positive 
constant so that ht > ^(1 + x'^)'^hxx- Let UN{x,t) = Ea; \Yt\ A A^, where the 
index indicates that Yq = x. By a maximum principle argument, compare 
[7], U]\fix,t) < h{x,t) independently of A^. Consequently, by monotone 
convergence, Eq \Yt\ < h{0,t) < oo. 

Theorem 12.31 thus applies and provides a generalised diffusion Xt which 
is a martingale started at such that Xi has the same distribution as Yi. 
Consequently, we have in this case two different generalised diffusions with 
the same starting point that give rise to the same distribution at time 1. 
However, only one of these processes is a martingale. 

6. An application to Mathematical Finance 

In this section we study an inverse problem in Mathematical Finance. Let 
X be a non-negative martingale with Xq = xq, and consider the expected 
values 

C{K, T) := ¥.{Xt - K)+. (6.1) 

Here Xt has the interpretation as a price process of a stock, and the expected 
value C{K,T) is the price of a call option with strike K and maturity T 
(for the sake of simplicity, we assume that interest rates are zero). The 
function C is non- increasing and convex as a function of K, and it satisfies 
C(0,r) = xo, C{oo,T) = and C{K,T) > {xq - T)+. 

In Mathematical Finance, the corresponding inverse problem is of great 
interest. Option prices C{K,T) are observable in the market, at least for a 
large collection of strikes K and maturities T, whereas the stock price model 
is not. Under some regularity conditions (often neglected in the literature), 
Dupire [5] determines a local volatility model 

dXt = aiXt,t)dWt 

such that (16. ip holds for all K and T. To do this, naturally one needs call 
option prices C{K,T) given for all strikes K >0 and all maturities T > 0. 

The assumption that call option prices are known for all maturities T is 
often unnatural in applications - indeed, there is typically only one maturity 
a month. We therefore consider the situation where C{K,T) is given for all 
strikes K but for one fixed maturity T > 0. A natural question is then 
whether there exists a time-homogeneous diffusion process 

dXt = a{Xt) dWt 

such that (16. ip holds for all strikes. We have the following result. 

Theorem 6.1. Let c : [0,oo) [0, cxo) he a convex and non-increasing 
function satisfying c(0) = xq, c(oo) = and c{K) > (xq—K)^, and letT > 
be fixed. Then there exists a time-homogeneous (generalised) diffusion such 
that X is a martingale with Xq- = xq, and such that c{K) = E(Xt — K)'^ . 
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Proof. We construct a probability distribution fi on [0, oo) as follows. On 
(0,00), we let fi be defined as the measure given by the second derivative 
of the convex function c. Since c{oo) = 0, this measure has mass — c'(0), 
where c'(0) is the (right) derivative of c at 0. Moreover, we let a point mass 
e = 1 + c'(0) be located at 0; note that the assumptions c(0) = xq and 
c{K) > {xo - K)+ >xo-K imply that c'(0) > -1 so e > 0. 

Integration by parts shows that the expected value of the distribution fi 
is given by 

x/x(dx) = / fj,{x,oo) dx = — / c'{x) dx = c{0) = Xo- 
/o Jo Jo 

Now, by Theorem 12.31 (and an obvious scaling to general T) there exists a 

time-homogeneous generalised diffusion X with Xq- = xq such that X is a 

martingale and Xj- has distribution fi. Finally, integration by parts gives 

E{Xt-K)+= {x-K)fi{dx)= n{x,oo) dx = c{K), 
Jk Jk 

thus finishing the proof. □ 
7. MARTINGALITY of GENERALISED DIFFUSIONS 

The speed measure constructed in the proof of Theorem 12.31 is such that 
the generalized diffusion (Xt) is a martingale, but as mentioned above, this 
is not the case for every speed measure v. We characterize here the speed 
measures for which Xt is a (local) martingale. 

Example 7.1. If the speed measure is given by 

dx for X > 



v{dx) 



for 2; < 0, 



then the corresponding process X is a Brownian motion reflected at 0, and, 
in particular, not a local martingale. 

This type of reflection at an extreme point of supp v is the only thing that 
can prevent X from being a local martingale. This is shown in the theorem 
below from here somewhat extended. 

We first give a lemma, essentially saying that stopping X at a point 
a € suppi^ is the same as stopping B aX, a. (Note that we cannot stop X at 
a ^ suppi^ since Xt G suppz^ for all t > by Lemma l3. 11 ) 

Lemma 7.2. If a & {x'^,x^) Hsuppz^ and r := mi{t : Xt = a}, then 
T = ■ If Ha < -ffsuppoo then also r < 00 and = Ha, while t = 00 if 
Ha > Hsupp^u- Furthermore, F{Ha < i^suppoov) > and P(r < to) > for 
any to > 0. 



Proof. If Ha < Hsuppaau, then Th^ < 00 by Lemma [331 and thus Lemma 
and the assumption a G suppi^ imply that Lf^^+e > Th^ for every e > 0. 
Thus, if t = THai then At = Ha and Xt = = = On the other 
hand, \it< Th^, then At < Ha by (|2.3p so Xt = BAt o-- Hence, r = Th^- 
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If Ha > Hsuppaou, then Lemma [331 yields Th^ = oo. Lemma [331 yields 
further, for any t < oo, At < Hsupp^u < Ha so Xt = BAt 7^ a. Hence, 

T = 00 = Th,. 

For the final statement, assume that a < xq, say. Note that the as- 
sumption implies xq ^ suppoo i^, so x^ < xq < x'^. Choose b G {xo,x'^). 
Then z/[a,6] < oo. Let e := to/via, b]. (The case z^[a, = is easy.) It is 
easy to see (e.g. using excursion theory) that with positive probability both 
supj. -L|^ < e and Ha < H^. Thus Bs € [a,b) for all s < Ha, and then 

= = fa L'lj^vidx) < ev[a, b] = to. 

The claim F{Ha < -ffsuppoo j^) > follows from this, using Lemma [3?5t but 
it is easier to see it directly since Hg^pp^^, = ^ 

Theorem 7.3. {Xt)t>o- is a local martingale if and only if either suppz/ = 
{xq} (then, trivially, Xt = xq for allt) or both the following conditions hold: 

(i) suppoo n (— oo, Xq] ^$ or suppi/ n (—CO, 6] 7^ for all b < xq, 

(ii) suppoo n [xq, oo) ^ 9 or supp z/ n [b, oo) ^ for all b > xq. 

Moreover, if {Xt)o-<t<to is a local martingale for some to > 0, then 
{Xt)t>o- is CL local martingale. 

Proof. If supp u = {xq}, then Xt = Xq a.s. for all t by Lemma 



Suppose now that (i) and (ii) hold. If suppoo ^ H (— oo, xq] ^ 0, let a~ := 
x°° = sup{x < Xq : X £ suppoo z^} for each n; otherwise let a~ be a sequence 
of points in supp z^H (— oo, xq) with a~ \ — oo as n — >• oo. Define a sequence 



On ^ ^0 similarly, using (ii) 

Let Hn ■= Hr- +x. For each n, (Bt/\H„)t>o is a bounded martingale. 
Thus, if s < t, then E{BAtAHn | -FaJ = Ba^ah^, so {BAtAHjt>o- is a 
martingale for the filtration (Qt) = iJ~At)- 

Set T„ := r/f„_|_ = inf{r„ : u > Hn}. Tn is a stopping time for the 
filtration (Gt)- Further, since Hn < Hn+i we have T„ < Tn+i. 

If Tn < oo, then T}{„ = < oo and it follows from Lemma 13.21 that a.s. 
then Tu > Tn for all u > Hn, so At„ = Hn, hence, Sa^ah^ = ^a^aAt^ = 
^AtATn ~ -^tAT„- On the other hand, it r„ = oo, then At < Hn for every t, 
and thus BAtAH^ = BAt = = ^Mr„- In any case, thus BAtAH^ = XtAT„ 
a.s., so {XtAT„)t>o- is a martingale. 

Finally, we verify that T„ — )■ oo a.s. as n — >■ oo. Recall that Bh„ G 
{a~,a^}; we consider for definiteness the case when Bh^ = so Hn = 
H^~. By construction, either a~ € suppoo and then Lemma 13.51 implies 
that Th = oo a.s., or else o~ — >■ — oo and then H„- — )• oo and thus 
Th _ — )• oo a.s. In both cases we obtain T„ = Th _+ — >• oo a.s. This 
completes the proof that {Xt)t>o~ is a local martingale. 



For the converse, we assume that suppi/ ^ {xq}. Assume that (i) fails 



(The case when (ii) fails is symmetric.) Let ^_ := inf{x : x € suppz^}. 



Since (i) fails, ^_ > — oo. We may assume that ^_ < xq, since otherwise 



Lemma 13.41 shows that Xt is not a local martingale at 0- 
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Let Ti = , the hitting time for Bu of Since (i) fails, [^_,xo] H 
suppoo = so G {x°^,x'^) and Lemma 17.21 shows that with positive 
probabiUty Ti = < Hgupp^u, i-e., Bu hits ^_ before it hits suppoo z^; 
denote this event by £. Lemma [7.21 further shows that if £ holds, then there 
exists s < cxD (viz. s = Tti) such that Xs = We want to show that in 
this case there also exists t > s with Xt 7^ 

If suppj^ = then ^_ < xq by our assumption suppv 7^ {xq}, but 

then Xt is not a local martingale at 0— by Lemma [3.4[ We may thus assume 
that there exists another point X2 7^ in suppz^. Let T2 be the first hitting 
time of X2 after Ti. We study two cases separately. 

If < 00, then Lemma 13.21 implies that a.s. Lrj+e > ^7^2 for all e > 0. 
In this case, if t = Ftj; then At = T2 and Xt = Bt2 = X2 ^ 

If = 00, then = infju : Tu = 00} is finite. If [/ is a neighbourhood 
of Ba^, we can find e > such that Bu £ U for u £ {Aoo — £, ^00 + e), and 
thus = L^^_^ for x ^ U. Since Taoc+s — ^Aoo-e = 00, it follows 

that I'iU) = 00. Hence Ba^^ ^ suppooi^, and in particular Ba^ 
t — )• 00, At — > Aoo and thus Xt = BAt Ba^ 7^ so in this case, Xt 7^ ^- 
for all large t. 

Note also that Xt > ^_ for all t > by Lemma 13.11 If {Xt)t>o- is a 
local martingale, then Y given hy Yt = Xt — is thus a non-negative local 
martingale, hence a supermartingale. Therefore zero is absorbing for Y, see 
e.g. ^14. Proposition (3.4)]. But this contradicts our result above which says 
that with positive probability Yt first hits zero and then takes a larger value. 
Hence X is not a local martingale. 

For the final statement, it suffices to show that ¥{Ss < io) > for some 
(5 > since the argument above then shows that {Xt)t<to is not a local 
martingale. With ri = inf{t > : Xt = = T^i and T2 = mi{t > ri : 
Xt^C-} it thus suffices to show ¥{t2 < to) > 0. 

By Lemma 17.21 P(ri < to/2) > 0. Let A := T2 — ti and condition on £. 
Then A < cx) a.s. and the strong Markov property of B implies that A is 
independent of ri and, moreover, that P(A >a + b\ A>a)= P(A > b) 
for any a,b > with P(A > a) > 0; thus (conditioned on f ) either A = 
a.s. or A has an exponential distribution; in both cases P(A < to/2) > 0. 
Hence, P(r2 < to) > which completes the proof. □ 

Before giving the corresponding characterization for martingales, we give 
some lemmas. 

Lemma 7.4. If At < H^^^pp^y, then VAt = t- Consequently, for all t > 0, 

TAt=t^ r^eupp^ . • 

Proof. If At < -f^suppoo u, then for small e > we hav e At + e < -f^suppoo f and 
thus ^At+e < cxo by Lemma [3?5t which by Remark [23 Hii) yields = t. 



Furthermore, in this case t = Fa^ < ^Hsupp^. so F^j = t = t A F^.^pp^,. 

Since At < -ffsuppoo f by Lemma 13. 5[ the only remaining case is At = 
i^suppoo^.. In this case F^.^pp^, = Fa^ < t so t A Th.^^^,^, = ^H^upp^. = 
Fa,. □ 
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Lemma 7.5. For every measurable function / > 0, a.s. for every u < 



/ f{Bs)dTs= / f{x)Llu[dx). (7.1) 

Jo J-oo 



Proof. By monotone convergence it suffices to consider u < Hgupp^u, and 
then, by modifying u outside the range of {Bg : s < u}, it suffices to consider 
locally finite v. In this case, the result is a consequence of the general 
theory of continuous additive functionals, see [141 Corollary X.(2.13)]. (It 
is also easy to make a direct proof in this case, by first assuming that / is 
continuous, and then partitioning the interval [0, u] into N subintervals, and 
for each subinterval estimating the change of the difference between the two 
sides in (17.111: we omit the details.) □ 



Lemma 7.6. If f > is a deterministic or random measurable function, 
then a.s. for every T < Hg^pp^ ^, 

r f{s)drs= f{At)dt. (7.2) 
Jo Jo 

Proof. Consider a fixed uj in our probability space. By a monotone class 
argument (or by seeing the two sides of (I7.2p as J f dfiL and J f djiR for 
two finite measures ^l-,1^r on [0, oo)), it suffices to prove (I7.2p when / is 
the indicator of an interval [0, u) for some n > 0. In this case f^ f{s) dTg = 
Jq'^^ dVs = Ttau and, by Remarks [27 



m 



f{At)dt= / l{t<Tu}dt = TThTu = VTAu- □ 
'0 Jo 

If {Xt)t>o- is a martingale, then EXt/x,- = xq for every X-stopping time 
r and every t > 0. This means that 'KBt = xq for the i?-stopping time 
T = At /\ At-. We would like to have E Bt = xq also for other i?-stopping 
times T < At, not necessarily obtained by stopping X. However, this is not 
always possible, as is seen by the following example. 

Example 7.7. If suppz^ = {xq}, then Xt = xq for all t > 0—, so Xt 
is trivially a martingale. However, if further z^{xo} < oo, for example if 
v = 5x0, then for any t > and a ^ xq, F{BAtAHa = a) = IP(^^a < At) > 0; 
since BAtAHa ^ {a,xo}, this implies ESajAH, / xq. 

The next lemma shows that Example 17.71 is the only counterexample. 
(The trivial example shows that the lemma is not as trivial as it might 
look.) 

Lemma 7.8. Suppose that supp v ^ {^o} and that {Xt)t<tQ is a martingale 
for some to > 0. Then, for every stopping time T < Atg and every real a, 

It will follow from Theorem 17.91 and its proof that, more generally, E Bt = 
Xq for any such T, so the Ha is not really needed but it simplifies the proof. 



22 



ERIK EKSTROM, DAVID HOBSON, SVANTE JANSON, AND JOHAN TYSK 



Proof. We suppose that a < xq. (The case a > xq is similar and a = xq is 
trivial.) 

Suppose first that o G suppz^. Then Lemma [3. 2 1 implies that a.s. THa+e > 
for every e > 0. Hence, if r := Tj^^, then = Ha and X-r = Bh^ = 
moreover, Xs ^ a for s < Ha- In other words, r is the X-stopping time 
inf{s : Xs = a}. The assumption that {Xt)t<to is a martingale thus implies, 
for t < Iq, 

xo = E Xt;,r = E BA,^Ar = IE ^A, a//. • (7.3) 
If T' and T" are two stopping times with T' > T", then ¥,(^Bj''/\u\J'T") = 

B 

T"Au for every u > (see e.g. [10, Theorem 7.29]). If T' < Ha, then 
— a > for all s < T', and Fatou's lemma yields K{Bj'' — a \ J-t") ^ 
Bxii — a and thus, by taking the expectation, 

¥.Bt'<^Bt". (7.4) 

We apply ([7^ first to At^ A Ha and T A Ha and then to T A Ha and 0, 
obtaining 

IESa,„ah„ < IE5tah„ < ESo = a;o, (7.5) 
and (j7.3p shows that we have equalities. This proves the result when a £ 
suppz^. 

In general, by Theorem l7.3l either there exists some b < a with b G suppz^, 
or there exists b E [a, xq] with b G suppoo i^- In the first case Ha < Hf, and 
in the second case Ha > Hi, and Lemma [3.51 implies that Af^ < Hi, and thus 
T < AtQ < Hj, < Ha- Consequently, in both cases T A Ha = T A Ha A Hi,, 
and the result follows from the case just proven applied to T AHa and b. □ 

Theorem 7.9. The following are equivalent, for any to > 0. 

(i) {Xt)t>o~ is a martingale. 

(ii) {Xt)o-<t<tQ is a martingale. 

(iii) suppi^ = {xq} or xq E suppoo or 

r-QO r-XQ 

j {l + \x\)v{dx)= j {l + \x\)v{dx) = oo. (7.6) 

J X(, J —oo 

(iv) supp V = {xq} or Xq E suppoo v or 

r-QO rxo 

/ \x — XQ\i'{dx) = / \x — xo\i'{dx) = oo. (7-7) 

J Xq J —oo 

Remark 7.10. For a related result for non-negative diffusion processes, see 

in- 

Proof. The result is trivial if suppz/ = {xq} or xq E suppoo z^- We may thus 
assume suppi^ ^ {^o} and — oo < x°^ < xq < x^ < oo. The equivalence of 
(|7.6p and ()7.7p then is elementary. 
Let 

{2 fr N \v\ y(dv), X > Xq, 
or / : (^-^^ 
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and 

ijj{x) : = 



n,^{y)dy, x>xo, 
fr'p{y)dy, x<xo. 



Then '0 is a non-negative convex function on {x'^,x'^) with left derivative 
(/? and thus second derivative (in distribution sense) ip"{x) = 2\x\v{dx). 
By the Ito-Tanaka formula \\.4\ Theorem VI. (1.5)] and Lemma 17.5^ for 

i^[Bu)= / ip{Bs)dBs+ / Ll\x\u{dx)= / ip{Bs)dB,+ / \Bs\dTs. 
Jo Jr Jo Jo 

(7.10) 

Let H := H^a,b} where x'i^ < a < xq < b < x"^. Then H < Hg^pp^^. 
Further, f{Bs) is bounded for s < H, and consequently J^^^ (p{Bs) dBg, 
u > 0, is a martingale. Hence, for every bounded stopping time T < H, 
E ip{Bs)dBs = 0, and fTTOD yields 

E(^^(St))=e/ \Bs\dTs. (7.11) 
Jo 

Moreover, for any stopping time T < H, we can apply (j7.1ip to T Au and let 
It — >■ oo. Since T < H we have Btau & [a, b], and thus iP{Btau) is uniformly 
bounded; hence dominated convergence on the left-hand side and monotone 
convergence on the right-hand side shows that (j7.1ip holds for any stopping 
time T < H. 

We apply (j7.1ip first to T = Aj. A H Au for some r,u > 0. Thus, using 



Lemma 17.61 and Ty^^ < r, see Remarks \2Z Ifiii) 



rArAHAu r-ArAHAu 
Hi'iBA^AHAu)) = E / \B,\drs=E \BsAHAu\ dV, 

Jo Jo 

r^ArAHAu rr rr 

= E \BAtAHAu\dt <E \BAtAHAu\dt= E\BAtAHAu\dt. 

Jo Jo Jo 

(7.12) 

Suppose first that ()7.7p holds. By translation we may assume that xq = 0. 
Then ip{x) ±oo as x ^ ±oo, and thus iIj{x)/\x\ og as x ^ ±oo. In 
particular, ip^x) > \x\ for large and thus |x| < Tp{x)+C for some constant 
C. Consequently, (j7T^ yields 

E \BArAHAu\ <C + E{iIj{Ba^ahau)) <C + [ K \Ba,ahau\ dt. 

Jo 

Note further that BAtAHAu is bounded by max{|a|,5}, so the expectations 
here are finite and bounded. We can thus apply Gronwall's Lemma [14^ 
Appendix §1] and conclude 

E\BA,AHAu\<CeK (7.13) 
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Using this in ()7.12p again yields the sharper estimate 

E(V(5a,aha«)) < CeK (7.14) 

Now let a and h x^; then H — )• ^^suppoo;^ = H{x°°,x'^}- Since 

At < Hsnpp^v by LemmaESl BAtAHAu BAtAu, and Fatou's lemma yields 

E(^(Ba,a«)) < Ce\ (7.15) 

Since 'il^{x)/\x\ ^ oo as x ^ ±00, (j7.15p implies that for a fixed t, the 
random variables BAtAm u >0, are uniformly integrable ^ Theorem 5.4.3]. 
Moreover, BAtAm u > 0, is a martingale, and since it is uniformly integrable, 
BaiAu = IE(i?At I J^u), for all u > 0, and further, see [lOl Theorem 7.29], for 
any s < t, 

E{Ba,\TaJ = Ba,aa. = Ba,. (7.16) 
Hence Xt = -Ba* is a martingale when (|7.7p holds, so 



(iv) 




(i) 


— ^ 


(ii) 



Conversely, suppose that (ii) holds but (17. 6p fails; by symmetry we may 
assume that 



{1 + \x\) iy{dx) < 00. (7.17) 

I XQ 

In particular, this shows that suppoo H (xo,oo) = 0, so x^ = 00. By 
translation, we may now also assume that x'^ < < xq, and we now take 
a = and any b > xq. Thus H = H^^ i^y = Hq A Hf,. We apply ()7.1ip with 
T = At A H. Noting that H < Hq and 2?^ > for s < H, we obtain, using 
also Lemma rr6t 

pAtAH i-AtAH 

E V'(Sa,aHoAhJ = E / BsdVs = E BsAH, dVs 

Jo Jo 

= E BA^AHodr. (7.18) 

Jo 

Since H < Hgupp^ v , Lemma 17.41 shows that 

Ta.ah = La, a Lh = t a Lh^^pp^. ArH = tArH = tArH,A Th,. 

Further, if r > Th^ then > Hq by (12. 3p and thus BArAHo = Bhq = 0. 
Hence, (17181) yields 

Eip{BAtAHoAHt) =^ BArAHodr = E BArAHo dr. 

Jo Jo 

(7.19) 

Let C := 2 y v{dy). Then |9j(x)| < C for x > 0, and thus, for any 
X > 0, 

V'(x) < V'(O) + C'x < ip{xo) + C'xo + C'x = C'xo + C'x. (7.20) 
Thus, for t = to, the left-hand side of (I7.19P is at most, using Lemma EHJ 
E(C"xo + C'Ba.aHoah,) = 2C'xo. (7.21) 
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Consequently (I7.19P yields 

2C'xo > E 



toAFj 



BArAHo dr. 



(7.22) 



Letting 6 — )• oo we have //^ — ?• oo and Thi, oo, and thus (|7.22p yields by 
monotone convergence, using Lemma 17.81 again. 

rto rto 
2C'xo > E / BArAHo dr = E BArAHo dr = IoXq. 
Jo Jo 



Consequently, 



This shows that 



to<2C' = A yu{dy). 



(i) 




(iii) 





(7.23) 



since we then may choose to arbitrarily large 



in (j7.23p and obtain a contradiction. To extend this to (ii) with a given to 



we argue as follows, assuming (ii) and (|7.17p 



We have derived (|7.23p under the assumption x'^ < < xo- By transla- 



tion, we have in general, for any xq and x"^ , and any a € {x^,xq) 

POO 

to < 4 / yi^idy). 

J a 

Letting a — ?■ xq yields 



to <4 



yv{dy). 



(7.24) 



(7.25) 



Xo 



Take any z > xq with z G suppi/, and let r := Th/, by Lemma [7.21 
T = inf{t : Xt = z}, so r is an X-stopping time. Condition on the event 8 := 
{t < to/2}; we have F{£) > by LemmaEH On the event £, B'^ = Bh,+u 
is a Brownian motion starting at z, and the processes corresponding to F, A 
and X defined by B' are = F^^^+^-F^^ = Th,+u-t, A[ = Af+r-H^ and 



Xi 



Xt+r- Since r is an X-stopping time, on X[ is a martingale 



for < t < to/2, and also at t = 0— since X'q = z = Xq_, cf. Lemma 
Consequently, we may apply the result above to X' , and (|7.25p yields 



to 
2 



< 4 



yv{dy). 



(7.26) 



However, if (j7.17p holds, then the right-hand side of (j7.26p tends to as 
z — > oo, which is a contradiction. This contradiction shows that (j7.17p 
cannot hold, and thus (ii) =^ (iii) □ 



Note that the case xq € suppoo v is not redundant in (iii) An example is 
given by Xo = and v{dx) = dx/x for x > 0; then (I7.6P fails because the 
second integral vanishes, nevertheless = xq = for all t, which trivially 
is a martingale. 
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8. Open problems 

8.1. Weakening the assumptions. Theorem 12.31 shows the existence of a 
generahsed diffusion solving the inverse problem provided the given distri- 
bution has finite expectation. Does the result hold without this assumption 
(allowing local martingales instead of only martingales)? 

8.2. Uniqueness. Theorem 12.31 shows the existence of a generalised diffu- 
sion, defined by a speed measure with a given distribution of Xi (as- 
suming E|Xi| is finite). Is the measure v unique? Strictly speaking it is 
not: Lemma 13.71 shows that we will never leave the interval so 
any change of v outside this interval (assuming at least one of x"^ and x"^ 
is finite) will not affect Xt at all. We may normalize v by first taking the 
restriction to {x'^,x'^); if x°^ is finite but x'^ ^ suppoo I'o we also add 
an infinite point mass at x!^', and similarly for xj^. The question is then: 
Given the distribution of Xi, is there a unique corresponding normalized 
speed measure such that the corresponding process {Xt)t<i is a martingale? 

The problem of uniqueness is also open in the discrete case, i.e. is the 
mapping G : — > H" studied in Section H] an injection? (As remarked in 
the proof, this holds for n = 1.) 

8.3. Relations between fi and z^. Find relations between the speed mea- 
sure (assumed normalized as above) and the distribution fi of Xi. For 
example, if n has a point mass at x, does also have a point mass there? 
Does the converse hold? If v is absolutely continuous, is fi too? Does the 
converse hold? For which i/ does fj, have a finite second moment? 
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